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Models of multi-field inflation exhibiting primordial non-Gaussianity have recently been intro- 
duced. This is the case in particular if the fluctuations of a light field scalar field, transverse to 
the inflaton direction, with quartic coupling can be transferred to the metric fluctuations. So far in 
those calculations only the ensemble statistical properties have been considered. We explore here 
how finite volume effects could affect those properties. We show that the expected non-Gaussian 
properties survive at a similar level when the finite volume effects are taken into account and also 
find that they can skew the metric distribution even though the ensemble distribution is symmetric. 



I. INTRODUCTION 



The observations of Cosmic Microwave Background anisotropics and of the large-scale structure of the universe offer 
a window to the physics of the inflaton. For instance the detailed measurement of the shape of the power spectrum 
can give constraints on the shape of the inflaton potential. On the other hand the detection of Non-Gaussian metric 
fluctuations could signal the existence of effective couplings between the inflaton and other fields. There have been 
a series of observational progresses in view of constraining the deviation from Gaussianity using CMB data on large 
scales (0 - 10°) Q, intermediate (0 ~ 1°) and small scales (0 ~ If/) 0, as well as large scales structures 
The recent CMB data by WMAP were also used by Komatsu et al. that concluded that the CMB spectrum 
was compatible with Gaussianity on the basis of an analysis of the bispectrum. Similar analysis using Minkowski 
functionals ||| and 3-point correlation function Q reached the same conclusions. More recently, an analysis of the 
WMAP data using measurements of the genus and its statistics 3 conluded that the Gaussianity of the CMB field 
was ruled out at a 99% level. As stressed after, non-Gaussianity of primordial origin may still have escaped detection 
and further investigations are needed. 

From a theoretical point of view, Gaussianity is a generic prediction of slow-roll single field inflation 0, TTol . A 
series of models including feature in the inflationary potential > existence of seeds such as x 2 A2\ , axion |l3j or 
topological defects [l4| , curvaton scenario 0| , varying inflaton decay rate ^(| have been shown to be able to generate 
some primordial non-Gaussianity. 

In Ref. jlCj . we proposed a general mechanism to produce non-Gaussianity in the adiabatic mode and explicit 
realizations of such models were presented in Ref. 01 • The mechanism is based on the generation of non-Gaussian 
isocurvature fluctuations which are then transferred to the adiabatic modes through a bend in the classical inflaton 
trajectory. Natural realizations were shown to involve quartic self-interaction terms. The statistical properties of the 
resulting metric fluctuations were then shown to be the superposition of a Gaussian and a non-Gaussian contribution 
of the same variance, the relative weight of the two contributions being related to the total bending of the trajectory in 
field space. The non-Gaussian probability distribution function (PDF) was also computed and shown to be described 
by a single new parameter so that generically only two new parameters suffice in describing this class of models. 

In order to infer constraints from e.g. CMB data on these kind of models predicting non-Gaussianity one needs at 
least to go through 2 steps. 

1. One first needs a precise prediction of the statistical properties of the curvature in order to construct estimator 
adapted to the detection of this kind of non-Gaussianity. For instance the analysis Komatsu et al. p| uses 
the bispectrum and Minkowski functional and constrains a x 2 deviation from Gaussianity. The gravitational 
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potential was parametrized as $ = <I>l + /nl( ( E > l — (^l)) where is the Gaussian linear perturbation of zero 
mean and it was concluded that —58 < /nl < 134 at 95% CL. This constrains only a very peculiar type of 
non-Gaussianity and does not apply e.g. to the models of Refs. [lfl Il7j. 

2. Even if the form of the PDF is predicted theoretically, one needs to investigate what is measured and how the 
measurements are related to the theoretical predictions. 

The goal of this article is to investigate the observational implications of the theoretical predictions of Ref. 01 ■ It 
requires further investigations in particular on the effect of finite volumes of the survey while the detailed comparison 
of our model to existing dataset, beyond the scope of this work, is still left for further studies. In particular, we want 
to show that these finite volume effects do not suppress the predicted non-Gaussianity. 

Section [H] formulates the problem, focusing on the quantities that can actually be observed. In Section IIHI we 
recall basics of the models we consider in this paper and give simple consequences. We then develop in Section IIVI 
some arguments to understand the effect of the mean value of the non-Gaussian component on the size of the survey. 
SectionEJis devoted to the computation of the lowest order measured cumulants. In particular, we will show that there 
appears a non vanishing skewness that is directly induced by the finite volume effects as stressed in our concluding 
remarks, Section IVll 



II. AN OVERVIEW OF FINITE VOLUME EFFECTS 



Basically what we want to measure are the statistical properties of a scalar field \ of zero mean 

<X> = 0, 

where and from now on angular brackets, ( ), refer specifically to ensemble averages. Such a scalar field will be in the 
following identified with the gravitational potential but the actual observations may be complex linear transforms 
of that field such as the temperature anisotropies or polarization of the CMB or the large-scale cosmic density 
fluctuations. The general questions raised by precision measurements of statistical properties of cosmic fields can 
be very intricate. We refer for instance to[lg Section VI for extensive developments regarding such issues. We will 
restrict our discussion to the overall aspects of them. 

In practice while performing such measurements, two scales enter the problem, the scale R s at which x is measured 
(e.g. somehow smoothed) and the size R n of the survey. The smoothed field Xs can be measured at different locations 
in the survey. The idea is that its spatial fluctuations should provide us with hints of the actual statistical properties 
of the field. Ensemble averages are however inaccessible as such. The values of x s we have access to are in finite 
number and are all inherited from the a single stochastic process, the one of our universe. One can however get 
insights into the stochastic properties of x from geometrical averages. In the following, we denote such a geometrical 
average A and its connected part A " . For instance one can have access to quantities such as 

x n = ( Xs - x) n ' = £x7 r °- 

The X n are themselves stochastic quantities that depend on the stochastic field x- To be more specific we have 

X=f x(y,i)W« H (y)d 3 y (1) 

and 

Xs W = J X(y, t)W Rs (x - y)d 3 y (2) 

where VVr is a window function of volume unity and for simplicity assumed to be spherically symmetric, e.g. for a 
top hat filtering we have VVr = 0(|x — y| — R)/V where G is the heavyside function and V the volume of the ball or 
Wr = exp(— |x — y| 2 /(2i?) 2 )/(-\/27ri?) 3 for a Gaussian window function. 

Clearly quantities such as (x 2 ) cannot be estimated because they contain contributions from modes of wavelength 
larger than the size of the survey and that cannot be observed. Only modes with wavelength in a given range are 
measurable. This is in particular the case for a Gaussian field. In this case one can show that if the two scales R s and 
i? H are in large enough ratio, X n can be good estimates of ((x s — x)"), e.g. the expected dispersion of the measured 
values X n decreases with i? H / R s to a power that depends on the power spectrum shape. In our case however not only 
the super-Hubble modes cannot be measured but the modes that can be observed are correlated and also correlated 
to the super-Hubble modes. 
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In the following, we will explore the consequences of both the fact that only modes in specific wavelength are 
measurable and that they are correlated. Ideally, cosmological models should be able to predict the PDFs of the 
measured cumulants, 

but such predictions are obviously difficult to do in general. We will see in the following how we can estimate these 
distributions for some families of multi-field inflationary models. 



III. MODELS 



A. Self-interacting scalar fields in de Sitter space 



Before we start to investi gate the statistics of the observed quantities let us recall the basics of the models we have 
in mind. In the models [T(l Il7l|. the non-Gaussianities are first generated by an auxiliary field, Xi self-interacting in 
a potential, typically quartic, 



V( X ) 



4! 



X 



(3) 



This field is assumed to be a test field so that it does not affect the dynamics of inflation driven by an other scalar 
field, <f>. Assuming an almost de Sitter inflation and neglecting the gravitational back-reaction on the evolution of the 
universe, its evolution will be dictated by 



a z 3! 



(4) 



where H — a/a is the Hubble constant and A the comoving Laplacian. The evolution of the scale factor is given by 

-Ht 



a(t) = aoc 



HI 



V = 



H 



(•5) 



where t is the cosmic time and rj, which is negative, the conformal time. 

For a free massless scalar field (A = 0), the Klein-Gordon equation Q in the de Sitter background takes the form 



v = 



where we have introduced v = x/ a - The general of this equation is 

Hrj 



Xk(t) 



>2k 



1 



1 

ikrj 



-ikr) 



once the quantum field, Xt is decomposed in Fourier modes as|2J 

d 3 k 



(2tt) 3 / 2 



Xk(t)e 



h.c. 



where [S k , SJ.,] = (27r) 3 <5(k - k'). 



It implies that the field x has a correlator given by 

( X (k)x*(k')) = (2^) 3 P(fc) ( 5(k-k') 



(6) 



(7) 



(8) 



(9) 



where P(k) is the power spectrum, which is equal to P(k) = H 2 /2k 3 on super-Hubble scales for a free scalar field 
living in de Sitter space. This is the so-called Harrison-Zel'dovich spectrum. 

Moreover, the self interaction term of x induces non-zero high order correlation functions. For a quartic potential 
the odd order correlation functions vanish; the even order ones can be computed from a perturbation theory approach 
at tree order. For instance, the four-point function in the super-Hubble limit reads 



<Xki ■ ■ -XkJc - 
as explicitly shown in Ref. j2|J . 



~(2nf l ° g{v ^ ki) S^K) P(h)P(k 2 )P(k 3 ) + perm. 



(10) 
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B. Second and fourth moment 



When one wants to relate the expectation values of observable quantities such as the X n to the statistical properties 
of the field, filtering effects should properly be taken into account. The observable quantities, 5x s , can be written in 
terms of the Fourier modes, Xki as 

6 Xs = J j0j5e tk x x(VW(k) (11) 



with 



It follows immediately that 



W(k) = W(kR s ) - W(kR H ). (12) 



(Sxl) = 4n J k 2 P{k)W 2 {k)Ak. (13) 

Such an expression can be easily computed in case of a Harrison- Zel'dovich spectrum and a fc-space top-hat window 
function. Defining fc s and fc H respectively as 1/R S and 1 / R H , the second moment of Sx s reduces to 



<^> = 2^ 2 m(£)^. 



(14) 



It follows from this calculation that the expectation value of the observable quantity X 2 is given by 

(X 2 ) = (5 X 2 s ) = al (15) 

Its (cosmic) variance, (Xf ) , can similarly be computed and it scales like R S /R H in case of a Gaussian field |29| . 
Furthermore, the expectation value of X4 is given by 



(X*) = (5xt) c = I — ; ;; 6 V (fco . . . w(k 4 ) ( Xkl . . . xk 4 ) c . (i6) 



d 3 ki . . . d 3 k 4 

Expressing the 4-point correlator by mean of Eq. I|1U|) and assuming a Harrison- Zeldovich spectrum and a top-hat 
window function in fc-space, this expression finally reads 

A\ f J-f^ \ P^s H 3 W f^s H 3 k" f^s H 3 k" 

(^>c=-^h- /. V/ V/ V ^ (|ki+k2+k3|)log[(fci+fc2+fc3+|ki+k2+k3|)??] - 

_ ^ (17) 

where W (|ki + k2 + k3 1), defined in Eq. (|12|l . simply expresses the condition k u < |ki +k2 +k3| < fc s and is inherited 
for the term W{k/C). When the ratio k s /k H is large this expression can be easily computed after noting that, in this 
case, one of the ki generically dominates over the others (the contributions to the integral being roughly uniform over 
the integration domain when the wave- vector norms are logarithmicallly spaced) . In this limit it is therefore possible 
to approximate ki -I- k2 + k3 by say ki and k\ + k 2 + k% + |ki + k2 + k3 1 by 2 k\ . Finally the integral reads, 

(^) c --^Alog(2^/ 4 fcy 4 ) ± (18) 

which corresponds to what has been found in our previous work provided the number of e-folds is identified with 
log (2?7 fc 3 / 4 fc^ 4 ) ■ It means that the kurtosis of the <5xs ne ld is significant if A log {2r\ k^ 4 k 1 / 4 ) approaches unity. 
This result demonstrates that the mechanism described in our previous study survives finite volume effects. We can 
also note that observable quantities are insensitive to the behavior of P{k) in the small k limit. 

Finite volume effects have however other consequences due to the fact that the observable modes are correlated. 
The physical reason of these correlations is that they share the same history, e.g. they have been produced in the 
same stochastic process. We found that a random walk approach for the evolution of the local x values gives 
precious insights on those more subtle finite volume effects. 
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IV. LESSONS FROM A RANDOM WALK APPROACH 



From the previous definitions, it is clear that x an d the different values of x s share contributions of super-hubble 
modes. Those cannot be observed but they shape the values of x an( I Xs ■ Actually the value of x R at a given scale, 
whether it is i? H or R s , is dynamically built from the stacking of modes that successively leave the horizon. 

A random walk approach can then be used to describe the stochastic growth of x R ■ It will allow two things, to get 
insights into the excursion values of x an d to see how x s values are correlated through their common history. In this 
approach the field value evolution is described in term of a Langevin eq during inflation. 

Before we go to this equation, let us sort out how the different scales and evolution equation and related together. 



A. Scales intervening in the problem 



Different scales will have to be distinguished in our study. 

• During the inflationary stage, the super-Hubble modes of the scalar field can be treated as classical. In a de 
Sitter spacetime, the physical Hubble radius is constant i?(P h y s ) = H^ 1 so that the comoving smoothing scale is 
time dependent, R = (aH)^ 1 . The evolution equation of the classical part will thus contain a stochastic force 
simulating the effect of the quantum noise due to the modes which are crossing the horizon at each time step to 
become classical. We thus define a stochastic field, x H M i f° r which the filtering scale is dependent and always 
equal to the horizon size, i?? 1 ^ 8 = a(t)R H with i?^ 8 - iJ" 1 so that 

X H = XR=l/aH- (19) 

The dynamics of x H will follow a Langevin equation. We investigate the dynamics and the statistical properties 
of x H m Section llVBI 

• On the other hand the field values x an d Xs can be seen as a time dependent quantity but they correspond to 
the filtering of x a t fixed physical scales. They identify though with x H a t precisely the time r\ at which the 
scales R H and R s respectively cross the horizon, e.g. x(v — — V-^ho) = Xh(v = — V-^ho) an< ^ XsC 7 ? = — V-^s) — 
X H (v = — V-^s)- After that these coincidental times the two fields x ancl Xs behave classically, i.e. follow an 
inflationary classical Klein-Gordon equation without stochastic source terms. 

To summarize, two of the comoving scales R s and R H0 are fixed while one, i? H is a time dependent quantity. A 
sketch of the different sequences that the field dynamics follow are shown on Fig. 2] 

B. The late time PDF of Xaiv)- 

We follow a formalism first developed to deal with self-interacting fields in a de Sitter background 0] based on 
the idea that the infrared part of the scalar field may be treated as a classical spacetime dependent stochastic field 
satisfying a Langevin equation [2(], 0] (see Ref. for the case of a massless free field) . 

Assuming that x 1S slow rolling, the dynamics of x H can be obtained by averaging Eq. (0J in which both the second 
time derivative and the Laplacian term can be neglected, since x H contains only long wavelength modes. Using the 
identity 

X H k = (xu) RH W(kR H ) - kR aX ^W'(kR u ) (20) 
where W'(u) = dW(u)/du and where we have used Eq. JSJ to express the time derivative of i? H . It follows that 



X H =-^'(x H +^x) + £ Q (x,t) ( 21 ) 



which reduces (5x ^ Xh) to 



*" = -^^ + ^'>- <22) 
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FIG. 1: The different filtered quantities and scales entering our problem. Xh follows a stochastic dynamics that can be 
described by a Langevin equation; it has a time dependent smoothing scale so that more and more modes contribute to the 
filtered field. The field values x an d Xs ev °l ve according a classical Klein-Gordon equation; their smoothing scale is time 
independent; they coincide with x H a ^ horizon crossing times. 



The term £ Q (x, t) appears from the commutation between (%) r h and x H using Eq. I|2U|) . It is a stochastic noise 
describing the effects of the small wavelength (quantum) part exiting the horizon on the classical stochastic part. 
Using Eq. ©, one obtains its expression as 

£ Q (x,t) =-J —^{kHR^) [xk(t)W'(kR H )e^b k + h.c.j . (23) 

Indeed, £ Q is quantum noise so that we replace it heuristically by a Gaussian stochastic noise, £ with a correlator that 
matches the quantum expectation value in the standard Bunch-Davies vacuum, that is 

(£(x, i)£(x', f ) = <0|£ q (x, i)£ Q (x', |0) . (24) 

It follows that 

- / ^MW' (^) W> , (25, 

with r = |x — x'|, which reduces, using the expression J7J on small scales, to |2^| 

(£(x,i)£(x',t')> = [ dfcsinfcr(l + ikr))(l - ikr}')^' -^W' (-kr))W' (-kr)') . (26) 

\ r n z r J 

Even though x H remains a quantum operator, it was replaced by a stochastic field in a way that, for all observables, 
the expectation values of the two fields are in excellent agreement. 

In the following and for the sake of simplicity, we consider a window function reducing to a top-hat in Fourier space 
so that W reduces to a Dirac distribution. In that case, using the solution (0) one obtains 

The case of more realistic window functions (such as a Gaussian) was discussed in Ref. [23|. 
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From a Langevin equation of the form 

X H =-pV' + a£, (Z( X> t)Z(x',t')) = 5(t-t') (28) 

where we assume the two coefficients and a to be independent of \ H > one can deduce |24j a equation, the Fokker- 
Planck equation, for the PDF of x H °f the form 

d t V=Pd* H V + ^d x jV'V]. (29) 

It follows from the Langevin equation (1221 with a top-hat window function for which the noise is given by Eq. Ij27(l 
with (3 = -1/3H and a = H 3 / 2 /2n that the one-point PDF, V{x H ,t), is solution of 

d ^-^lV + ^d x JVV). (30) 

In a cosmological context, this equation was first derived in Ref. |l9l | in the case H — const., which we are interested 
in, and then in Ref. j2^] in the case V = 0. It was generalized to more involved situations in Ref. plf . 

The solution of Eq. (|30|l where studied in Ref. [21| in which it is shown that V approaches the static equilibrium 
solution 

7> oq = AT 1 exp (-g^( X „)) , M = |" exp (-Uf(x h )) d XH (31) 

irrespectively of the initial conditions. 

For the quartic potential, Eq. |J3}, we find 

v M = 2rwm{~) exp 

They are few lessons to learn from this result. The excursion values of x H are bounded. Their distribution does 
not depend on the remote past history of x and, importantly for the following, the typical value one can expect for 
X H , and therefore x, is H/X 1 / 4 . 



7T A 

9iJ 4 



(32) 



C. Consequences on the shape of the PDF of 8x s 

We can then insights into the shape of the probability distribution function of <5x s as a function of x- 
Quantitative results can be drawn from the perturbation approach we initially developed in our previous work [lOfl . 
We expand the filtered field in terms of the coupling constant as 

xM=X ( °Hr,)+x£Hv) + ---- (33) 

X^ -* represents the value of the filtered field when the interaction term is switched off. In the slow-roll regime the field 
X s follows the Klein-Gordon equation Contrary of the previously studied case of x H ; t ne evolution equation does 
not contain any noise term because the smoothing scale R is now fixed and there is no new modes entering x s ■ The 
free filtered field, X$i is therefore constant and, from the discussion of the previous section, it cannot be assumed to 

be Gaussian distributed with a zero mean: its expectation value, Xs > is actually the value of x a t time r\ — l/i? H . 
Since this field value is going to be only weakly affected but its subsequent evolution, in the following wc will identify 

Xs '' an d X- The difference ^x^ = X^ ~ xf ^ nas been built up from the modes that have left the horizon between 
— l/fc H and — l/fc s . It can be assumed to be Gaussian distributed with a width precisely given by a$. 
The first order term in A, x s ^\ evolves according to 

3ffXs (1) =~[xi 0) ] 3 . (34) 

In this approach the treatment of the filtering of the r.h.s. of this equation is very crude. The results we are going 
to find can anyway be checked against more rigorous calculations based on the computed shape of the tri-spectrum. 
Our goal now is to capture the essential effects of a non- vanishing x on the statistical properties of 6x s ■ 
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FIG. 2: PDF of S\ s = 5\ s — x f° r different values of X- The dashed line corresponds to a Gaussian distribution; the dot-dashed 
line to the deformed distribution of 5x s when \N e /H 2 = 1 and X = an d the solid lines when the latter equals 0.5 and 1. 



The equation of evolution l|34(l can be solved to get 



(xf) 3 

18 H ' 



which also reads 



XN e 
18 H 2 



N e being the number of e-fold between tu and the end of inflation. 
These results imply that 



X ~ X 



(0) 



XN e 
18 H 2 



(0) 



o (0) 2 

3 Xs 



which explicitly shows that x and Xs°^ are equal at leading order in A. It also gives 



Sx s = <5xl 0) 



18 H 2 



(^ 0) ) 3 + 3[(<5xi 0) 



<y% ]x + 35xi 0) x 2 



It is straightforward to see that Sx s has acquired a non zero third order moment, 



XN, 



-X° '8 



(35) 



(36) 



(37) 



(38) 



(39) 



at leading order in A. This is a finite volume effect in the sense that it exists for a fixed (not ensemble averaged) value 
of x- This effect cannot a priori be neglected. From the study of the previous paragraph, we know that x should be 
of the order of -ff/A 1 / 4 , which implies that the reduced skewness of <5x s , fixl/i^xl) 3 ^ 2 ~ X 3 ^ 4 N e ag/H is significant as 
soon as X 3 ^ 4 N e approaches unity, a condition similar to that encountered in Sect. IIII Bl 



Actually the evolution equation for x s can be solved in 



Xs = 



x( Q) 



1 



A N , 
9H 



a ( (o)Y 

f (Xs ) 



(40) 



and the distribution of y, can then be infered from that of xi assuming the latter is Gaussian distributed with a 
non-zero mean value |30| . 

In figure we present the deformation of the PDF of <5x s while x is varied. As expected, it shows that when x is 
not zero, the PDF gets skewed in a way that can be easily understood: when x is positive it gets more difficult to 
have excursion towards larger value of x s > but easier to roll down to smaller values. It is as if the field x w &s actually 
evolving in the potential A(x + x)/4!- As a result one naturally expects the field x to have non-vanishing three-point 
function. As mentioned before, these calculations treat the smoothing in a rather simplified way but we think that 
for illustrative purposes it encapsulates the main effects that we want to describe. 

It nonetheless shows the way for the computation of the finite volume effects for expected stochastic properties of 
the field. 
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V. FINITE VOLUME EFFECTS ON HIGH ORDER CUMULANTS 



What the Langevin picture suggests is that finite volume effects on the observed quantities are not due to the whole 
stochastic process that give birth the observed field but mainly to the value of x alone. In other words, we expect to 
have 

T(X n ) = J V(x) S (X n - (( Xs - XT) X ) d*. (41) 

This form l|41(l implies in particular that the PDF of X n , V(X„), is expected to be peaked around the expectation 
value of 5x~ at x fixed. In particular, it implies that 

ra-((<^s),n ( 42 ) 

We will explicitly check such a property for the lower order cumulants. 

In general, for small enough values of x, the constrained ensemble averages of the form (A(x))^ should be given by 

{A( X ))- ~ (A( X )x) 4 (43) 
a x 

where cr| is the variance square of the x fluctuations. This relation, exact for Gaussian fields, is only approximate in 
general. It can be derived for stochastic variables following a quasi-Gaussian distribution. Here, it will be valid only 
if the excursion values of x are modest compared to the fluctuations of $x s • 
Not surprisingly, it implies that the even order cumulants are left unchanged. 



A. The bi-spectrum 



Non trivial finite volume effects are then going to appear at the level of the third order cumulants or correlation 
functions. In particular it induces a non-vanishing bispectrum, three-point correlation function of the wave-vectors, 
which is going to be given by 



( X (k 1 )x(k 2 )x(k 3 )) x * (x(k 1 )x(k 2 )x(k 3 )x) c 4 



(44) 



when x is small enough. Using Eq. (fTT)f> to express (xk a • ■ ■ Xk 4 ) c and using that Eq. J3J) implies that Xk = XkW(kRy), 
the previous expression reduces to 



(X(k 1 ) x (k 2 )x(k 3 )) x = -^(2^) 9/2 J, log 



p(fc 1 )p(fc 2 )p(fc 3 )+p 



E k < 



E k * 

i 

] [P(fci)P(fc 2 )+perm.] \W 



E k * 



p H 



(45) 



The factor W (|ki + . . . k 3 |P H ), that arises from the contribution of modes with k < 1/P H to x, ensures that |ki + 
. . . k 3 1 is small compared to each of the fcj , and thus can be neglected in the log term. It implies that for an Harrison- 
Zel'dovich type spectrum that P(|ki + . . . + k 3 |) ^> P(h L ) so that the first term of l|45|l is negligible. As a result we 
deduce that 



(x(k 1 )x(k 2 )x(k 3 )) i 



Ax 
3P 2 



(2^) 9 / 2 log 



J2 k 



Mb 1 P (IE^|)W/(IE^|P H ) (AM 
[P(fci)P(fc 2 ) + perm.J 5 • (46) 
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Now, noting that, by definition, P(k)W(kR H )d 3 'k/a^ integrates to unity and that the function 

P (Ei W (Ej ^1 ^h) i s J f° r the modes we are interested in, peaked near the origin we obtain that this 
factor is essentially equal to 5 (£\ kj). It therefore implies that 



( X (k 1 )x(k 2 )x(k 3 )) i 



Ax 
3H 2 



(2tt) 9 / 2 log 



[P(fci)P(fc 2 )+perm.] <H E k * 



(47) 
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Here is one of the main points of this paper: Finite volume effects induce a non- vanishing three-point function although 
the potential in which \ evolves is symmetric. 

From this bispectrum it is possible to compute the third order moment of 6\ s ■ Its amplitude will be in agreement 
with what had been obtained from the Langevin equation. It is also the three-point function one expects for a field 
evolving in potential A%x/3! (see Ref. [2(j for details). 



B. The three-point cumulant 



We now turn to the lowest order cumulant exhibiting a non-trivial result due to a non- vanishing x, that is X3, 
which reads 

X 3 =J d ^ )9 /2 k3 ^(l k i +k 2 +k 3 lfl H )^(fc 1 )^(fc 2 )^(fc 3 ) X k 1 yk 2 Xk 3 - (48) 
Obviously its ensemble average vanishes, 

(*3> = 

but not (Xf) . Let us check, as expected from our analysis [see Eqs. l|4H42fl ]. that it is well approximated by 

To evaluate the latter expression, we start from the expression of (Xa)^ that is defined by 

W s *§sj ^^g/f^ ^dki + • • ■ + k 3 \R H )W(h) . . . W(k 3 ) < Xkl • • . x* 3 )t ■ (49) 
Using Eq. Q47jl. it reduces after integration over IC3 to 

(X 3 ) t ~ J d 3 k 1 d 3 k 2 M?(fc 1 )TU(A :2 )W?(|k 1 + k a |) {P(fci)P(fc 2 ) + P(|kj + k a |) [P(fci) + P(fc 2 )]} (50) 

where the window function W(|ki + ... + k3|P H ) was aborbed during the integration over k3 due to the Dirac 
distribution. Such an expression can be computed following the same lines as for the computation of the expectation 
value of the fourth order cumulant l|17fl . 

Ay f H 2 \ 2 r k s d 3 ki f k s d 3 k, ~ 
(X3> ^"3^ T / V/ -T^^dki+^D log [(fei + fe 2 + 1^+^1)7?]. (51) 

where, again, W (|ki + k 2 |) simply expresses the conditions, k H < |ki +k 2 | < fc s . A simple expression for this integral 
can be obtained when k s is much larger than k u where it is possible to replace ki + k 2 and k\ + A; 2 + |ki + k 2 | by 
respectively either ki and 2 k\ or k 2 and 2 fc 2 . Finally the integral reads, 



<*3> x - -Ax(4tt) 2 f £L) log 2 ( j log fakgX") (52) 



e.g. 

T 4 



(X 3 )^ = -X X log (2^ 2 / 3 fcV 3 ) ^ (53) 
which reproduces the result (|39|l if N e is identified with log (2?yfc 2/ ' 3 /c ] l i //3 ) . In conclusion, we end up with 

(^3),) 2 ) 1/2 -Alog(2^„ /3 )# (54) 



In this case, it is actually possible to compute .Xf ) from a perturbation theory approach in order to check that its 
dominant contribution is indeed < ( (X 3 ) - ) V (Xfy is given in general by 



t?(|kx + . . . + k 3 \R H )W(\k[ + ... + k' 3 \R H ) ( Xkl . . . Xk 3 Xk< 1 • ■ • Xk^) c • (55) 
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It involves the expression of the 6-point correlation function (xki • • ■ Xk-jX^ ■ ■ ■ Xkf,) c which, in a perturbation theory 
approach, can be split into 2 contributions 



T _ I (o) (o) (l) (o) (o) (l) 



where 



This implies that 



and 



£ W(x) 



18 P 2 



^ (0) (x) 



(2) (x) 



2 Ar2 r 



8 P 4 



C (°J(x) 



/ = ( 27 r) 9 I ^18 2 ^|P(fc 1 )P(fc 2 )P(k3 + k x + k 2 |)P(fc 2 )P(fc 3 ) 6 (J2 k * + E k 



1/ = (2^) 9 - 6! ^P( fcl )P(fc2)P(fci)P(fc 2 )P(fc^) 5 (V k, + £ k; 



(56) 



(57) 



where 18 2 and 6! are symmetry factors. 
Let us evaluate the first contribution 



(X 2 )^ = A 2 ^ / d 3 k 1 d 3 k 2 P(fc 1 )P(fc 2 )W r (fc 1 )^(fc 2 ) / ^kaPdka+kt+kaD'WC&a^dki+ka + kali^) 



2 JV e / j3, i3 



P 4 

d 3 k;d 3 k , 2 P(fc 1 )P(fc^)VK(fc' 1 )^(fc^)W ? (|k 1 + k 2 + k 3 + k[ + k 2 | 



(58) 



The term W 2 (\ki + k 2 + k 3 |P H ) implies that |ki + k 2 + k 3 | < P H so that W(k 3 ) ~ Wflki + k 2 |) and Wflki + k 2 
k 3 + ki + k' 2 | ~ W(\k[ + k' 2 |). Setting e = ki + k 2 + k 3 , we conclude that 



P 4 



d 3 k 1 d 3 k 2 P(fc 1 )P(fc 2 )H/(fc 1 )VF(fc 2 )W^(|k 1 + k 2 |) 



d 3 eP(e)l¥ 2 (eP H ). 



The integral over e reduces to <r| so that 



2\« 



N 2 

\2 e 2 

P 4 X 



d 6 k 1 d 6 k 2 P(k 1 )P(k 2 )W(k 1 )W(k 2 )W{\k 1 + k 2 |) 



(59) 



(60) 



The second contribution reduces to 



N 2 

90\ 2 j^ I d 3 k 1 d 3 k 2 P(fc 1 )P(fc 2 )^(fc 1 )V^(fc 2 )^(|k 1 +k 2 |) 



d 3 k' 1 d 3 k 2 d 3 k 3 P(fc 1 )P(A ;2 )P(A ; 3)^(fc 1 )Ty(A :2 )iy(fc 3 )^ 2 (|k 1 + k 2 + k 3 |P H ). 
The second integral reduces to 

d 3 k' 1 d 3 k 2 d 3 P(^)P(fc 2 )P(|k' 1 + k' 2 \)W(k[)W(k' 2 )W(\k' 1 +k 2 |) / d 3 eiy 2 (eP H ). 



TV 2 / - - 

90\ 2 -^R- 3 I d 3 k 1 d 3 k 2 P(fc 1 )P(fc 2 )V^(fc 1 )iy(fc 2 )^(|k 1 +k 2 |) 



(61) 



The integral over e gives ~ P H 3 so that 

r d 3 k' 1 d 3 k 2 P(fci)P(fc 2 )P(|k' 1 + k' 2 \)W(k[)W(k' 2 )W(\k' 1 + k' 2 \) 



(62) 



Due to the term W^(|k^ + k 2 |), we deduce that, in the case of a Harrisson-Zel'dovich spectrum, (R S /R H ) 3 < \k[ 
k 2 |~ 3 P^ 3 < 1 so that this contribution is at most equal to the one of Eq. ffity . As a result we have 



2\C0 



(63) 
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From Eq. I|60|l . this reduces to 

(Xl) 1/2 c\N e (2n)^ (64) 

that can be identified with the expectation value of (<5Xs)^ over the distribution of \, as obtained in Eq. J53J. 

This explicit computation shows that, as expected, the fluctuations of the measured values of X3 are mainly due 
the fluctuations of x- R justifies for instance that one should expect to see a bispectrum of the form (|T7|l for such 
inflationary models. 



VI. CONCLUSIONS 



In this article we have focused on the phenomenology of the non-Gaussianity generated in models developed in 
Refs. [ToL Il7| . Interestingly, whereas the metric perturbation statistics involve only two microscopic parameters 
respectively related to the weight of the non-Gaussian component and to its PDF, the finite volume effects imply 
that the statistical properties of any observational quantity will involve a third parameter. This new parameter arises 
from the fact that the mean value, x, of the field on the size of the observable universe does not vanish a priori. 
Obviously \ cannot be determined on the basis of any observation. As described in ^3 II VI it implies that the originally 
symmetric PDF can be skewed and that this skewness is directly proportional to \ that needs to be considered as a 
new parameter of the PDF while dealing with observations. 

These results open the way for more detailed phcnomenological studies. To see how those properties translate for 
the temperature and local density fields is not easy. Such investigation will probably require numerical tools such 
that developed in Ref. [27) . 

Acknowledgements: We thank Nabila Aghanim and Tristan Brunier for discussions as well as Alain Riazuelo for 
early discussions on the project. 
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